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SUMMARY 


/ 

J-j 

In  considering  life  tesC  experiments  from  a  behavioristic  Bayes 
point  of  view,  we  may  speculate  on  the  effect  of  model  mlsspeciflcation 
on  the  posterior  mean  under  alternative  stopping  rules.  Of  course, 
under  the  life  distribution  model,  we  expect  the  mean  of  the  posterior 
to  equal  the  mean  of  the  prior.  However,  If  we  think  the  model  may 
have  been  mlsspeclf led,  we  would  no  longer  expect  this  equality.  Ue 
show  that  for  Bayes  estimators  of  mean  life  under  the  exponential  model 
and  general  priors  on  mean  life,  we  expect,  based  on  a  preposterior 
analysis,  that  the  mean  of  the  posterior  will  increase  with  sample  size 
when  the  true  model  has  an  increasing  failure  rate  and  certain  fixed 


stopping  rules  are  used.  Also,  we  show  that  for  Bayes  estimators  of 
mean  life  based  on  the  natural  conjugate  prior  for  the  exponential 
model  and  complete  samples,  the  Bayes  risk  is  less  when  the  coeffi¬ 
cient  of  variation  is  less  than  for  the  exponential  model.  Recommen¬ 
dations  relative  to  use  of  life  test  sampling  plans  based  on  an 


exponential  model  are  made. 


PREPOSTERIOR  ANALYSIS  OF  BAYES  ESTIMATORS  OF  MEAN  LIFE,  I 


by 

Richard  E.  Barlow  and  Alex  Wu 


1.  INTRODUCTION 


In  Chls  paper  we  consider  selected  li£e  test  sampling  plans  with 
respect  to  life  distribution  model  robustness.  Our  objective  Is  to 
obtain  results  useful  in  selecting  life  test  sampling  plans.  To  this 
end,  let  F(*  |  9  )  be  a  life  distribution  with  mean  life  0  e  0  ■ 

[0,»)  .  Since  8  is  unknown  to  us,  we  model  our  uncertainty  about  8 
by  letting  8  denote  a  random  variable  with  density  ir(0)  .  We  suppose 
ir  is  chosen  and  fixed.  It  is  important  to  remark  that  0  and  ir 
may  be  considered  independently  of  our  life  distribution  model  F(*  |  0) 
subject  only  to  measurability  restrictions.  In  some  parts  of  this  paper 
we  will  concentrate  on  the  prior  density 


ir<e) 


bV(a+1)e"b/0 

T(a) 


(1.1) 


where  a  ,  h  >  0  .  Later,  we  also  assume  a  >  2  .  This  is  the  inverted 
garma  prior  which  also  happens  to  be  the  natural  conjugate  prior  when 

F(x  |  0)  -  1  -  exp  (-x/0)  . 

With  this  setup,  suppose  we  choose  and  perform  a  life  test  experi¬ 
ment,  E  ,  and  observe  data,  D  .  After  specifying  a  life  distribution, 
F(*  |  0)  ,  and  prior  ir  ,  we  may  compute,  ir(0  |  D)  ,  the  posterior 
density  for  0  given  the  data,  D  .  Let  expected  mean  square  error 


be  our  criterion  for  merit  of  en  eetiaetor  of  8  .  Then  we  seek  8(D) 
such  that 


E  _{0  -  8(D)>2 

"if 

is  minimum.  Expectation  is  with  respect  to  the  probability  measure  on 
tha  sample  space,  determined  by  our  sampling  plan,  and  9  ,  Induced  by 
our  life  distribution  model  F(*  |  8)  and  prior  density  ir  .  It  is 
well  known  that 

8(D)  ■  E(9  |  D)  •  1 8ir(6  |  D)d8  ,  (1.2) 

the  mean  of  the  posterior  density,  ir(8  |  D)  ,  minimizes  expected  mean 
square  error  among  estimators  based  on  the  data.  E(8  |  D)  is  called 
the  Bayee  estimator  of  8  with  respect  to  squared  error  loss.  Also, 
it  is  well  known  that 


8(D)  •  E(0  |  D) 


is  Bayesian  unbiased;  i.e.. 


E  _{9(D)} 
™ » t 


E 


if.F 


( 


0*(8  |  D)d0 


8 

o 


(1.3) 


where  0O  Is  the  mean  of  the  prior  density,  it  .  That  is,  we  expect, 
preposterior  analysis,  the  posterior  mean  to  be  the  same  as  the  mean 
of  the  prior.  Since  we  have  not  actually  observed  any  data,  this  is 
perfectly  consistent.  Hence,  given  ir  ,  F( •  |  8)  and  our  experiment 
E  ,  the  mean  of  the  posterior  is  Bayesian  unbiased  and  has  least  expected 
mean  square  error  among  estimators  of  0  based  on  the  data,  D  . 
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Suppose,  however,  chat  we  are  committed  co  using  a  certain  life 
distribution  model;  namely  the  exponential 

G(x  |  9)  -  1  -  exp  (-x/9)  . 

This  often  happens  in  industry.  A  contract  may  specify  the  life  dis¬ 
tribution  model  and  procedures  to  be  used  to  demonstrate  a  reliability 
goal.  If  this  is  the  case,  we  may  speculate,  prior  to  observing  data, 
about  the  effect  of  alternative  life  distribution  models  on  our  esti¬ 
mator.  For  experiment  B  ,  let 

9g(D)  -  E(9  |  D) 

be  our  Bayes  estimator  for  0  where  the  posterior  is  computed  using 
the  exponential  life  distribution  model,  G(*  |  9)  .  If  G  is  in¬ 
correct  and  F(«  |  0)  is  a  more  accurate  model,  then  in  general, 

Vf'Vd)>  *  \  ■ 

i.e.,  our  estimator,  based  on  G  ,  is  not  necessarily  unbiased  if  obser¬ 
vations  come  from  F(*  |  0)  .  Also, 

ew  F(e  -  qf(d)}2  <  E^ye  -  £g(d)}2  d.4) 

so  that  neither  is  9„(D)  best  with  respect  to  expected  mean  square 
error.  Our  problem  is  to  evaluate  Che  extent  of  these  discrepancies 
with  respect  to  the  sampling  plan  chosen,  i.e.,  the  number  tested  and 
the  stopping  rule  employed,  as  well  as  for  alternative  classes  of  life 
distributions.  We  call  this  paper  a  preposterior  analysis  because  we 


are  in  effect  speculating  about  the  posterior  distribution  before  any 
data  are  actually  observed.  We  hope  to  be  able  to  recommend  which 
sampling  plan  to  employ  when  the  model,  say  the  exponential,  is  actually 
misspecified. 

At  the  data  analysis  stage  with  the  data  in  hand,  we  might  pursue 
a  Bayesian  robustness  study  relative  to  our  life  distribution  model. 

For  example,  we  might  consider  a  class  of  life  distributions  F(*  |  o, 0) 
indexed  by  mean  life,  9  ,  and  a  parameter  a  .  For  a  *  1  ,  we  might 
have 

F(x  |  o  -  1,9)  -  1  -  e"(x/6)  . 

Then,  for  a  joint  prior  ir(a,9)  ,  we  compute  the  marginal  posterior 
density 

ir*(9  |  D)  -J*ir(a,9  |  D)da 

and 

ir(9  |  D)  “  it(8  |  a  ■  1,D)  . 

If  the  posterior  densities,  or  their  means,  are  "close"  we  then  con¬ 
clude  that  our  procedures  based  on  a  ■  1  ,  and  our  data  D  ,  are 
robust  for  the  family  F(*  |  a, 9)  ,  Box  and  Tiao  (1972,  Chapter  3).  We 
do  not  pursue  this  type  of  robustness  here.  In  contrast,  our  specula¬ 
tions  and  analysis  are  all  in  the  mind.  Our  specified  prior,  it  , 
is  our  only  reference  with  respect  to  knowledge  about  the  mean  life, 


Our  main  results  concern  Bayesian  bias  of  6(D)  under  selected 
stopping  rules  when  we  conjecture  that  F(*  |  9)  may  actually  have 
decreasing  mean  residual  life  conditional  on  9  .  See  Theorems  3.1  and 
4.1.  Our  second  criterion  is  the  expected  mean  square  error  of  our 
estimator  when  F(*  |  9)  is  not  exponential.  Results  for  this  cri¬ 
terion  are  given  in  Sections  2  and  3.  This  criterion  is  explored  in 
greater  depth  in  Part  II. 


2.  COMPLETE  SAMPLES 


In  this  section,  we  suppose  n  units  are  life  tested  until  fail¬ 
ure  so  that  our  data  D  *  (x, ,x~,  x  )  consists  of  n  observed 

l  /  n 

lifetimes.  Given  F(*  |  9)  ,  these  are  assumed  independent  conditional 
on  9  .  If  the  random  lifetime,  X  ,  has  distribution 

G(x  |  9)  »  1  -  exp  (-x/9)  , 

1  n 

then  x  *  —  7  x.  is  clearly  sufficient  for  0  .  Let  it  be  the  natural 
n  “  i 

conjugate  prior  given  by  (1.1)  with  a  >  2  .  Then,  the  Bayes  estimator 
is 

9(D)  ■  (1  -  w)9  +  wx  (2.1 

o 


where 


9q  -  b/(a  -  1)  ,  w  -  y2/{±  E  Var(X  |  0)  +  y2} 


and 


Y2  -  b2/ (a  -  l)2 (a  -  2)  . 

2 

We  assume  henceforth,  that  a  >  2  so  that  y  <  *  .  For  the  prior 
given  by  (1.1),  w  ■  n/ (a  +  n  -  1)  . 

It  can  be  shown,  (Biihlmann,  1970)  that  (2.1)  is  actually  the  least 
squares  linear,  in  observations,  approximation  to  the  mean  of  the  pos¬ 
terior  density  for  any  F(*  |  9)  and  any  prior  tt  with  mean  8q  and 
2 

variance  y  .  For  the  exponential  distribution  G(*  |  0)  ,  (2,1)  is 


the  mean  of  the  posterior  only  if  ir  is  the  natural  conjugate  prior 
(Diaconis  and  Ylvisaker,  1979). 

Clearly,  for  any  F(*  |  9)  and  it  , 

E  _{9(D)}  -  (1  -  w)9  +  w9  -  9„  (2.; 

ir,F  o  o  o 

so  that  (2.1)  is  always  Bayes  unbiased.  Also,  for  any  F(*  |  9)  and 
prior  ir  , 

2 

E  _{9  -  9(D)}2  =  (1  -  w)2y2  +  —  E  Var (X  |  9)  (2.: 

ir,F  n 

where  X  has  distribution  F(4  j  9)  .  Since 

2 

lin  w  =■  lim  y - ^ - - - -  -  1  , 

n-*«  n-*«  —  e  Var(X  |  9)  +  y 

it  follows  that  the  limiting  expected  mean  square  error  is  zero  for  all 
F(*  |  9)  and  ir  .  Hence,  sample  size  n  can  be  chosen  to  achieve  any 
specified  expected  mean  square  error  given  F(*  |  9)  and  u  . 

Let  CV(X  |  9)  -  /Var(X|  9 ) / 9  ,  the  coefficient  of  variation  of 
X  .  Let  Y  have  distribution  G(x  |  9)  =*  1  -  exp  (-x/9)  . 

THEOREM  2.1.  If  CV(X  |  9)  £  (£)  1  for  all  9  e  0  ,  then 

Et  _{9  -  9(D)}2  <  (>)  Ew  p{9  -  9(D)}2  ,  (2.< 

ir ,  F  ir  ,G 

where  9(D)  is  given  by  (2.1). 

Proof.  For  the  exponential  distribution  G(*  |  9)  ,  CV(Y  |  9)  -  1  for 
all  9  .  Hence,  if  CV(X  |  9)  £  (_>)  1  ,  then 


8 


4. 


I 

} 


t 

V 

* 

t : 

'  t  v»  '■»  •**  ■ 


Var  (X  |  9)  <_  (>)  9  2 

implies 

E  Var(X  |  9)  <.  (>)  E92  -  E  Var(Y  |  9)  .  (2.5) 

(2.4)  now  follows  from  (2.5)  and  (2.3). || 

From  Theorem  2.1,  it  is  clear  that  if  F(*  |  9)  is  a  member  of 
any  of  the  classes  IFR,  IFSA,  NBU  or  NBUE,  see  Barlow  and  Proschan 
(1975),  then  the  expected  mean  square  error  for  9(D)  is  less  than  that 
for  the  exponential  model.  Hence,  the  sample  size  n  can  be  determined 
assuming  model  exponentiality  to  meet  an  expected  mean  square  error 
criterion  and  it  follows  that  the  expected  mean  square  error  will  be 
even  less  if  F(-  |  9)  is  in  any  of  the  previous  classes.  Of  course, 
if  F  is  DFR,  DFRA,  NWU  or  NWUE,  the  expected  mean  square  error  will  be 
greater  than  in  the  exponential  case. 

The  Bayes  risk  is  the  expected  mean  square  error  with  respect  to 
F(*  |  9)  and  w  when  9(D)  is  the  posterior  mean  based  on  using 
G(*  |  9)  .  From  (2.4),  we  see  easily  that  CV(X  |  9)  1  implies  the 

Bayes  risk  under  F(*  |  9)  and  ir  is  less  than  the  Bayes  risk  under  G 
since 

min  Ew  (6  -  a)2  -  Ew  {8  -  E(9  |  D)}2 
a  if.F  ir,F 

<  Ew  _{9  -  9(D)}2 
—  ir,G 

where  E(9  |  D)  is  computed  with  respect  to  the  posterior  distribution 
based  on  F(*  |  9)  . 


TH 

3.  OBSERVATION  UNTIL  THE  K  FAILURE 

In  chis  section,  we  suppose  that  n  units  are  life  tested  until 
the  observed  failure  occurs,  k  is  fixed  in  advance.  We  call  this 
sampling  plan  (a).  At  this  point,  testing  stops.  Although  this  sampling 
plan  is  perhaps  not  as  useful  as  the  plan  which  stops  at  a  preassgined 
time  t  ,  it  is  easier  to  study.  Also,  results  for  this  stopping  rule 
suggest  similar  results  for  the  case  of  time  truncated  sampling, 
sampling  plan  (b) . 

Our  data  consists  of  the  k  ordered  lifetimes  D  *  (x^  £  x^ 

***  <  x,,.}  .  If  the  random  lifetime  X  has  distribution 
-  (k) 

G(x  |  0)  -  1  -  exp  (-x/0)  , 


then  the  total  time  on  test,  T  *  £  x^^  +  (n  -  k)x^  ,  is  sufficient 

for  9  .  Let  ir  be  the  natural  conjugate  (1.1).  Then,  the  Bayes 
estimator  is 


0(k)  -  (1  -  w)0q  +  w  |  (3.1) 

2  m  „  2  j 

where  w  *  y  /<  —  E  Var(X  |  9)  +  y  r  .  Since,  under  exponentiality  and 

conditional  on  0  ,  T  is  the  sum  of  k  independent  exponentials,  (3.1) 

is  actually  the  least  squares  linear,  in  T  ,  approximation  to  the  mean 

of  the  posterior  density  for  any  prior  rr  with  mean  0q  and  variance 
2 

Y  as  in  Section  2. 

Although  (3.1)  is  Bayes  unbiased  for  the  exponential  model  and  any 
prior,  it  is  not  necessarily  unbiased  for  arbitrary  F(*  |  0)  . 
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THEOREM  3.1.  If  F(*  |  9)  has  decreasing  mean  residual  life  condi¬ 
tional  on  9  ,  then 

E_  -C0Ck)>  >  6ft  •  (3.2) 

ir,F  o 

The  prior  is  the  natural  conjugate  (1.1).  Compare  this  with  Theorem  4.1 
in  the  next  section.  The  inequality  is  reversed  if  F(*  |  6)  has  in¬ 
creasing  mean  residual  life. 

Proof .  Assume  F(«  |  8)  has  decreasing  mean  residual  life.  For 
k  »  n  , 


0(n) 


(a  -  1) 
a  +  n  -  1 


9 

o 


a  +  n 


and  _.9(n)  ■  9  .  Applying  Bayes  Theorem  and  properties  of  the 

TT  9  r  O 

natural  conjugate,  we  have  for  k  <  n  , 


9(n) 


(1  -  w)8(k)  +  w 


(3.3) 


A  Q  —  k 

where  w  «  — — - -r  and  T{x„  *,*•}  is  the  total  time  on  test  from 

a  +  n  —  i  (K; 

x^^  to  the  last  ordered  failure  time.  Taking  expectations,  we  have 


9 

o 


/a  ±  k  - 
\a  +  n  - 


t) 


E9(k)  + 


(rSf^r) 


E8*(x(k)} 


(3.4) 


*. 

where  9  { x 
survival  to 


is  the  mean  residual  life  conditional  on 
x(fc)  .  By  assumption,  8  lx^k)}  £9  for  all 


8  and 
X(k) 


given 


Er,F9*lx(k)}  ^9o 


and  9  so  that 
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Hence,  from  (3.4),  we  must  have 

E  *9  00  1  6  • 

if ,  r  —  o 

The  proof  when  F(*  |  9)  has  decreasing  mean  residual  life  is 
similar. | | 

We  call  the  class  of  failure  distributions  with  decreasing  mean 
residual  life  DMRL.  If  F(‘  j  0)  is  IFR,  then  F(*  |  0)  is  DMRL. 

This  is  not  true  in  general  if  IFR  is  replaced  by  IFRA.  However,  a 
useful  comparison  is  available  under  star  ordering  (<p  of  distribution 
functions.  See  Barlow  and  Proschan  (1975,  Chapter  4.) 

Remark:  If  F  <  G  and  G  is  exponential,  then  F  is  IFRA. 

THEOREM  3.2.  If  F^*  |  0)  $  F2(‘  |  0)  V0  e  0  ,  then 

K  v  <90*)}  >  Ew  (0(k)}  .  (3.5) 

0(k)  is  given  by  (3.1)  and  is  calculated  under  G  .  This  strengthens 
(3.2)  relative  to  star  ordering. 

Proof.  From  (3.1),  0(k)  ■  (1  -  w)0q  +  w  £  and 

w  -  y2/  E  Var (x  j  0)  +  y2} 


IC  follow*  from  Barlow  and  Proschan  (1966)  that 

E_  (T  |  0)  >  E_  (T  |  0)  . 

F1  F2 

Hanca,  E  _  (T)  >  E  _  (T)  .  Applying  this  rasulc  to  (3.1),  tha 
~  n,t2 

thaoram  follows. | | 

Under  IFR  assumptions,  Bayesian  bias  actually  increases  with  n 
for  fixed  k  .  We  prove  this  and  more  in 

THEOREM  3.3.  If  F( •  |  0)  is  IFR,  then 

E  p{0(k)  |  n} 
ir,r 

increases  in  n  for  fixed  k  and  arbitrary  prior.  0(k)  ■  E(0  |  k,T) 

2 

and  the  prior  is  arbitrary  with  mean  0Q  and  variance  y  .  Note  that 
8(k)  is  the  posterior  mean  under  G  and  arbitrary  prior.  It  is  not 
necessarily  linear  in  T  . 

Proof.  Assume  F  is  IFR.  It  can  be  shown  that  when  ir  is  arbitrary, 
the  posterior  mean  under  G 

E(0  |  k,T) 

is  increasing  in  T  for  fixed  k  .  In  Barlow  and  Proschan  (1966),  it 
is  shown  that  when  F(«  [  0)  is  IFR, 

P(T  >  t  |  k,n,0) 

is  increasing  in  n  for  fixed  k  and  conditional  on  0  .  It  follows 
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E  _{9(k)  I  #)  •  I  _(E(9  I  k,T,n) } 

"»■  " «  * 

la  incraaalng  in  n  . | | 

If  F(*  |  9)  la  OFt,  than  E.  _{9(k)  |  a)  dacraaaaa  In  n  for 

W  f  F 

flxad  k  .  For  Cha  natural  conjugata  prior  and  9(k)  glvan  by  (3.1), 

E  _{9(n) }  -  9  It  followa  that  If  n.  <  n,  and  F(*  |  9)  la  in, 
»,r  o  i  i 

than 

Ew  r{9(k)  |  n2>  >  Ew  F{S(k)  |  nj  >  E^yStk)  |  k  -  n)  -  9o  . 

All  thla  la  trua  for  flxad  k  and  tha  natural  conjugata  prior  (1.1). 


If  F(*  |  9)  la  DF»,  all  lnaqualltla 


3.1  Tha  Expactad  Maan  Squara  Error  Crltarlon 

For  aampling  plan  (a),  not  only  tha  Bayaaian  blaa  but  alao  tha 
axpactad  aaan  aquara  arror  can  ba  unfavorabla  for  9(k)  ,  glvan  by  (3.1). 
In  particular.  If  F(x  |  9)  haa  ordar  of  contact  at  tha  origin,  vlth 
raapact  to  x  ,  graatar  than  that  of  an  axponantlal  distribution 
G(x  |  9)  ,  than 


11a  E  _{ 9  -  9(k) }  •  - 

_ "  *  * 


for  flxad  k  .  For  axaatpla,  thla  la  trua  If 
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i  ; 


F(x  6) 


1  -  exp 


(-<* -  VI 
(  <e  -  VJ 


X  <  X 


*ixo  • 


(3.6) 


where  9  >  x  >  0  . 

—  O 


THEOREM  3.4.  Assuae  either  F(x#  |  9)  ■  0  for  xq  >  0  or 


11*  1  9)  -  c  >  0 


(3.7) 


x+0  x 


where  a  >  1  .  Thet  ie,  F(x  |  9)  Is  0(x“)  as  x  *  0  .  Then,  for 
9(k)  ,  given  by  (3.1)  and  fixed  k  , 


lln  E  _{9  -  9(k)}2  -  •  . 
*«* 


(3.8) 


Proof.  It  is  easy  to  verify  thet 


E  _{9  -  9(k)}2  -  Var  -  9(k)}  +  E2  _(9  -  9(k)}  .  (3.9) 

ir,F  x,r  "  >  * 


Hence,  it  is  sufficient  to  show  for  fixed  k  that 


li«  _{9  -  9(k)>  -  -  . 

-  "  t* 


Since  T  -  nx^j  +  (n  -  l)(x^j  “  X(j^))  +  •.•■*'  (°  ~  k  + 1)  (x^j  -  x^_^j)  , 
it  is  enough  to  show 


Suppose  F(xq  |  0)  •  0  sod  xq  >  0  .  Then,  obviously 

11b  nE(X,  |  8)  >  11a  nx  -  •  . 
in  “*  o 

Froa  condition  (3.7),  it  follows  that  F(*  |  8)  is  In  the  doaain 
of  attraction  of  a  Welbull  distribution  with  shape  parameter  a  >  1  . 
See  Barlow  and  Proschan  (1975,  p.  241).  It  follows  that 

lia  p|(cn)^aX.  >  x  |  al  ■  exp{-(Xx)°}  for  all  x  0 

it-H»  ‘  n  ' 

where  8  ■  r(l  +  l/a)/X  .  Hence, 


11a  E 


rd  +  i/o)/x 


8 


since  the  convergence  is  uniform  in  x  .  Therefore, 


lia  «E(Xln  |  8)  «  lia  n1"1/a8/c1/° 


since  a  >  1  by  assumption. | | 


For  general  F  ,  we  have 


(3.10) 


2v  Cov 


-  8<*»2  ‘  «2[v«,{'(*  I  •)}  *  I  •)}] 
tov,}'(l  I  i)-i\ *  v“,(i>  *  *2[‘.  E(t  1  *)  ‘  9oJ  ' 


Mow  let  F(x  |  8)  -  1  -  exp(-(Xx)a)  for  a  >  0  *nd  8  -  ^ t 


Fro*  (3.10)  and  for  k  -  1  ,  we  have  when  0  <  a  <  1  , 


lim  E  _{8  -  8(1) }2  -  Var  (6)  +  wV  . 

It  ,1  IT  O 


IC  follows  that  for  k  ■  1  and  the  Weibull  distribution  with  o  >  0  , 


11*  E  _(8  -  8(1)}2  >  li*  E  _{0  -  0(1) )2  . 

r~-  “n—  *’C 


Hence,  asymptotically,  for  k  •  l  ,  the  expected  mean  square  error  for 
the  Weibull  distribution  is  always  greater  than  or  equal,  asymptotically 
to  the  expected  mean  square  error  in  the  exponential  case. 

Suppose  k  •  n  ■  1  ,  then  for  o  >  1  , 


E  _{0  -  0 (1) )2  <  E  _{0  -  8(1) )2  . 


Hence,  for  a  >  1  end  k  -  1  ,  the  expected  mean  square  error  as  a 
function  of  n  is  initially  less  than  the  corresponding  expected  mean 
square  error  for  the  exponential  case.  It  can  be  shown  that  eventually 
for  some  n  it  is  greater  and  remains  greater  thereafter.  In  the 
limit,  of  course,  it  is  infinite.  Theorem  3.4. 


4.  OBSERVATION  UNTIL  TIME  t 


In  this  section,  we  suppose  n  units  are  life  tested  until  time 
t  .  This  is  a  common  and  important  life  test  sampling  plan,  called 
sampling  plan  (b).  If  the  random  lifetime  X  has  distribution 

G(x  |  6)  «  1  -  exp  (-x/0) 

k 

then  k  ,  the  observed  number  of  failures  and  T  ■  J  x,n  +  (n  -  k)t  , 

i-1 

the  total  time  on  test  are  together  sufficient  statistics  for  9  . 

For  the  natural  conjugate  prior  (1.1),  with  a  >  1  ,  the  mean  of  the 
posterior  is  now 


0(t)  - 


(1  -  w)9  +  w  t-  for  k  >  1 

ox  — 


(4.1) 


9  +  — ~ 

o  a  -  1 


for  k  ■  0 


where  w  ■  k/(a  +  k  -  1)  .  Under  exponentiality  and  the  natural  con¬ 
jugate  prior. 


E  r{0(t)}  -  9 
ir,G  o 


(4.2) 


by  (1.3)  since  the  posterior  mean  is  always  Bayesian  unbiased  when  ex¬ 
pectation  is  computed  with  respect  to  the  model.  However,  in  general, 

E  _,{ 0  (t ) }  +  9  • 

it,  r  o 

THEOREM  4.1.  If  F(*  |  9)  has  decreasing  mean  residual  life,  i.e., 
is  DMRL,  then 


( 
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t 


E.  _{0(t) }  18  Vt  >  0  . 
ir,r  o 


(4.3) 


it  Is  the  natural  conjugate  (1.1)  for  the  exponential.  The  Inequality 
Is  reversed  If  F(*  |  9)  has  increasing  mean  residual  life. 

Proof.  Assume  F(*  |  9)  has  decreasing  mean  residual  life.  If 
t  ■  •  ,  then  we  have  complete  observations  and  our  estimator  becomes 

0  ■  (1  -  w)0  +  wx 

'  o 


and  E  _9  ■  9  for  all  F(*  I  9)  in  this  case,  see  Section  2.  Now 
ir,r  o 

suppose  k  failures  occur  in  [0,t)  .  By  Bayes'  Theorem  and  proper¬ 
ties  of  the  natural  conjugate  prior. 


<l  -  i>S(0  ♦ 


9(t) 


k  <  n 


(4.4) 


where  w  »  — ^ - — =■  and  T(t,<")  is  the  total  residual  life  of  the 

remaining  n  -  k  units  at  time  t  . 

Let  *  (k,T(t)}  be  the  observed  data  for  the  interval  [0,t)  ; 
l.e.,  the  number  of  failures,  k  ,  and  the  total  time  on  test  T(t)  in 
[0,t)  .  Since  F(*  |  9)  has  decreasing  mean  residual  life  conditional 
on  9  ,  we  have 


E  Jw  j  D. }  -  E  { |  D.i  <  ,  ^  ~  r  E(0  |  D  )  .  (4.5) 

ir,Fl  n-k  1  1)  ir,F(a  +  n-l  1  1)  —  (a  +  n-1)  1  1 


Hence, 


e„,f{“  I  Dl}  i  SE<5  1  Di>  ■  ”®<c>  • 


Note  that,  given  D1  ,  the  distribution  of  9  is  always  updated 
using  G  .  Hence,  E(0  |  -  9(t)  .  Also,  w  -  0  if  k  -  n  . 

It  follows  that 

9  -  E  _8  <  E  r{(l  -  w)9(t)>  +  Ew  -,{w9(t)} 

o  n,F  —  ir,F  ir,F 


and 

0  1  E  9(t) 

o  —  ir,r 

as  was  to  be  shown. 

All  inequalities  are  reversed  when  F(*  |  0)  has  increasing 
mean  residual  life. | | 

Remarks .  If  F(*  |  9)  is  IFR,  then  it  is  also  DMRL  so  that  (4.3) 
holds  in  this  case.  Since  F(*  |  9)  IFRA  does  not  in  general  imply 
F(*  |  9)  is  DMRL,  we  do  not  know  at  this  time  whether  or  not  F(*  |  9) 
IFRA  implies  (4.3). 

The  classical  sample  theory  approach  to  the  exponential  model 
and  sampling  plan  (b)  is  very  unsatisfactory  (Bartholomew  1957,  1963, 
and  Barlow  and  Proschan  1967).  The  principle  difficulty  is  that  if  no 
failures  occur  in  tO,t)  ,  the  MLE  is  infinite.  The  sample  distribu¬ 
tion  properties  of  the  MLE  are  also  highly  unsatisfactory.  In  the 
Bayesian  approach,  the  Bayes  estimator  (4.1)  is  Bayesian  unbiased  under 
exponent iality.  There  is  no  sample  theory  analogue  to  Theorem  4.1. 
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5.  CONCLUDING  REMARKS 


For  complete  sampling,  there  Is  no  Bayesian  bias  using  (2.1). 

Also,  the  expected  mean  square  error  Is  finite.  For  distributions  in 
the  classes  IFRA  and/or  NBUE,  DFRA  and/or  NWUE,  the  preposterior 
expected  mean  square  error  is  less  than  under  exponent iality.  Hence, 
for  complete  sampling  and  F(*  |  8)  IFRA  and/or  NBUE,  we  may  proceed 
under  the  exponential  model  assumption.  We  can,  by  appropriate  choice 
of  sample  size,  in  this  case  control  the  preposterior  expected  mean 
square  error. 

However,  for  the  sampling  plan  which  stops  at  the  kth  observation, 
or  time  t  ,  we  have  positive  Bayesian  bias  if  F(*  |  9)  is  DMRL. 

Hence,  if  it  is  necessary  to  use  stopping  rule  (a)  or  (b) ,  and  we 
have  strong  reason  to  believe  that  F(*  |  8)  is  not  exponential, 
then  we  should  consider  a  Bayesian  robustness  study  at  the  data  analysis 
stage,  see  Introduction.  The  more  IFRA  or  DFRA  in  the  sense  of  star 
ordering,  see  Barlow  and  Proschan  (1975,  Chapter  4),  the  more  compelling 
the  reason  for  a  robustness  study  at  the  data  analysis  stage  if  we  are 
committed  to  using  an  exponential  model  based  estimator. 
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